Non-Levy Distribution of Commodity Price Fluctuations 
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Price fluctuations of commodities like cotton and wheat are thought to display probability dis- 
tributions of returns that follow a Levy stable distribution. Recent analysis of stocks and foreign 
exchange markets show that the probability distributions are not Levy stable, a plausible result since 
commodity markets have quite different features than stock markets. We analyze daily returns of 
29 commodities over typically 20 years and find that the distributions of returns decay as power 
laws with exponents a which have values a > 2, outside the Levy-stable domain. We also find that 
the amplitudes of the returns display long-range time correlations, like stocks, while the returns 
themselves are uncorrelated for time lags m 2 days, much larger than for stocks ( ~ 4 min). 
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The study of economic markets has recently become an 
area of active research for physicists, for various reasons. 
First, markets constitute complex systems for which the 
variables characterizing the state of the system — i.e., the 
price of the goods, the number of trades, the number of 
agents are easily quantified. Second, because of the im- 
portance of markets there is a large amount of data that 
can be accessed. 

Much of the research interest of physicists has concen- 
trated on stocks JlJ , stock averages (5) , foreign exchange 
rates || and other markets Q . A number of key empir- 
ical findings have been established: (i) The distribution 
of logarithmic price changes is approximately symmet- 
ric and decays with power law tails with an exponent 
a + 1 « 4 for the probability density function 
(ii) The price changes are uncorrelated beyond rather 
short time scales and (hi) the amplitude of the 
price changes have long-range correlations, specifically, 
the correlations decay as a power law with an exponent 
7 « 1/3 0. 

Perhaps the most intriguing aspect of these empirical 
findings is that they appear to be universal, that is indi- 
vidual US stocks appear to conform to these "laws" Q , as 
do German stocks [||, and Australian stocks @. More- 
over, market indices such as the S&P 500, the Dow Jones, 
the NIKKEI, the Hang Seng or the Milan index || also 
obey these same laws. Furthermore, similar results are 
found for the most traded currency exchange rates such 
as the US dollar vs. the Deutsch mark, or the US dollar 
vs. the Japanese yen ||. The "universal" nature of the 
statistics of daily returns is remarkable since the mar- 
kets described above are quite different in their details. 
The observed universality thus suggests a rather striking 
similarity in the underlying mechanisms. 

Unlike stock and foreign exchange markets, commod- 
ity markets have received less attention fy-]lj]. Con- 
trary to heavily traded stocks or currencies — which have 
a somewhat abstract character because they (i) have an 
almost "elastic" response to changes in demand, (ii) do 
not require storage, and (hi) are not "consumed" — com- 



modities are physical products that are traded because 
they (i) cannot be produced at will, (ii) require physi- 
cal storage, and (hi) are needed for some purpose. For 
example, one needs gasoline to run a car, heating oil to 
heat a home, or electricity to light an office. 

Since stock markets and commodity markets differ in 
so many respects, one might imagine that commodity 
prices show larger fluctuations than stock prices. In fact, 
probability distributions of the returns of commodities 
such as cotton and wheat have been reported ||-|ll| to 
be Levy stable that have power law decaying tails with 
exponents in the range < a < 2, which implies the 
probability of a large gain or loss is larger than if the 
distribution were a power law with a > 2 0, || . Here we 
address the question whether the scaling of commodities 
is statistically distinguishable from that of stocks. We 
study the statistical properties of the price fluctuations 
for 29 commodities [|15[ and compare the results with the 
statistical properties of daily returns in stock markets. 

We define the normalized price fluctuation ("return") 

where At = 1 day, i = 1,2, ..,29 indexes the commodi- 
ties, Si(t) is the price and Oi is the standard deviation. 
Figure |l] displays the price and corresponding returns of 
a typical commodity, high-sulphur fuel oil [Tabic fif . Fig- 
ure H shows that the probability distributions P{gi > x) 
for both positive and negative tails decay as power laws 
with exponents that we cannot distinguish statistically, 

P(9i > x) ~ (2) 

where oti is outside the Levy stable domain < on < 2. 

Figure || displays Hill estimates Jl6| of on for all 
29 commodities, calculated for x above a cutoff value 
^cutoff) which is estimated by evaluating a, starting from 
different trial values of cc cu toff, an d choosing x C utoff as 
the minimum x value above which on does not change 
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significantly. Based on our analysis we choose for all 
29 commodities the same value x cuto ff = 2. Note that 
for daily returns the number of data points beyond 
a^cutoff = 2 is typically 50-200. The average exponent 
is a — (Si=i a i)/29 = 2.9 ±0.08 for the positive tail and 
a = 2.7 ± 0.07 for the negative tail 0. 

Next we compare our calculations of on with exponents 
on of daily returns evaluated for 7128 stocks from the 
CRSP database JT|,|l^]. We choose stocks in the same 
time period as the 29 commodities analyzed, and com- 
pute tail exponents of P(x) by the same procedure. Fig- 
ures ||(a) and |](b) compare the probability density func- 
tions of tail exponents for both commodities and stocks. 
The pdf for stocks is same as is reported in , and the 
exponents are outside the Levy stable region for both 
stocks and commodities. 

We next discuss time correlations of returns. Fig- 
ure ^(a) displays the autocorrelation function < 
gi(t)gi(t + t) > averaged over the 29 commodities. We 
observe that this averaged autocorrelation function C (r) 
ceases to be statistically different from zero for time lags 
r of 3 days or more. To further quantify time correla- 
tions, we use the detrcnded fluctuation analysis (DFA) 
method ]l9| ]. The DFA method calculates fluctuations 
F(n) in a time window of size n, and then plots F(n) 
versus n. The slope a in a log-log plot gives informa- 
tion about the correlations present. If C(r) ~ r~ 7 then 
a = (2 - 7 )/2, while if C*(r) - e" T / r - then a = 1/2 f§. 
We find in Fig. ||(b) that a — 0.55 ±0.05, consistent with 
the exponential decay of Fig. |S|(a). We also observe that 
|<7i|, the absolute value of returns (one measure of volatil- 
ity), are power law correlated with a — 0.73±0.05, which 
implies a power law decay of the autocorrelation of the 
absolute value of returns with 7 = 0.54 ± 0.1. 
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For stock markets it is seen that the autocorrelation func- 
tion decays exponentially to the noise level within a time 
scale of 4 minutes [Q,|^|. 
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TABLE I. Description of commodities analyzed 



Index 


Symbol 


Description 


Period 


# points 


1 


Brent 


Crude oil (UK Standard cf WTI in the US) 


1/88-8/98 


2770 


2 


BUTANE 


Butane chemical used as fuel 


2/93-8/98 


1433 


3 


Gasoil 


Gas oil 


1/88-7/93 


1433 


4 


HFO 


Heavy fuel oil 


1/88-8/98 


2770 


5 


HSFO arg 


High-sulfur fuel oil from Arabian Gulf 


1/88-8/98 


2770 


6 


HSFO New High-sulfur fuel oil transported by New York barge 


1/88-8/98 


2770 


7 


Kero New 


Kerosene transported by New York barge 


1/88-8/98 


2770 


8 


LSFO New 


Low-sulfur fuel oil transported by New York barge 


1/88-8/98 


2770 


9 


LSFO NYH Low-sulfur fuel oil traded at New York Harbour 


1/88-8/98 


2770 


10 


Nap Med 


Naphtha from Mediterranean (used for feed-stock) 


1/88-8/98 


2770 


11 


Nap New 


Naphtha transported by New York barge 


1/88-4/95 


1897 


12 


Prem unl 


Premium unleaded automobile gasoline 


6/92-8/98 


1619 


13 


CL 


Crude Oil 


1/83-8/99 


4164 


14 


HO 


Heating Oil 


9/79-8/99 


4999 


15 


C 


Corn 


6/69-8/99 


7593 


16 


CT 


Cotton 


6/69-8/99 


7593 


17 


FC 


Feeder Cattle 


7/79-8/99 


5051 


18 


KC 


Coffee 


6/69-8/99 


7593 


19 





Oats 


6/69-8/99 


7593 


20 


s 


Soybeans 


6/69-8/99 


7593 


21 


SB 


Sugar 


1/80-8/99 


4909 


22 


w 


Wheat 


1/75-8/99 


6186 


23 


LH 


Live Hogs 


6/69-8/99 


7593 


24 


PB 


Pork Bellies 


6/69-8/99 


7593 


25 


GC 


Gold 


6/69-8/99 


7593 


26 


HG 


Copper High Grade 


1/71-8/99 


7195 


27 


PA 


Palladium 


6/87-8/99 


3039 


28 


PL 


Platinum 


6/87-8/99 


3039 


29 


SI 


Silver 


6/69-8/99 


7593 
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FIG. 3. (a) Exponents on of the negative tail and (b) the 
positive tail, where i = 1, 2, .., 29 indexes the 29 commodities 
analyzed [Table §. We employ Hill's method |u| to estimate 
the exponent at of each probability distribution in the range 
x > Xcutofi where we choose x cn toH = 2. The dashed lines 
show the average values a = (X^ 2 =i a i)/29 . Shaded regions 
indicate the range of Levy stable exponents, a < 2. 



FIG. 1. (a) Prices Si(t) for a typical commodity, 
high-sulphur fuel oil HSFO arg [Table ^. (b) Normalized re- 
turns Qi(t) defined by Eq. (^). Note the large fluctuations 
of up to 20cr, which would have probability ~ e~ 200 for a 
Gaussian distribution. 
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FIG. 4. Probability density function of (a) the negative 
tail exponents and (b) the positive tail exponents evaluated 
for the 7128 stocks and the 29 commodities analyzed. Observe 
that for both stocks and commodities the mean exponent is 
around 3.0, outside the Levy stable region < a < 2. 



Normalized return, x 

FIG. 2. Cumulative distributions of positive and nega- 
tive returns for high-sulphur fuel oil (HSFO arg) [Table jj. 
The tails of the distributions shown have power law decays 
with an exponent as — 2.6 ± 0.3 for the negative tail and 
as = 2.9 ± 0.06 for the positive tail, where the exponents and 
the error bars are estimated by Hill's method Jl6| , 
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FIG. 5. (a) Autocorrelation function averaged over all 29 
commodities. The autocorrelation function displays an expo- 
nential decay with a decay constant of 2.3 days (see inset), 
(b) The fluctuations F(n) of returns <?i and absolute value of 
returns \gi\ averaged over all 29 commodities. The plot for 
gi is consistent with the exponential decay of Fig. 5(a). The 
exponent of a — 0.73 ± 0.05 for \gi\ implies long-range power 
law correlation, a feature also seen for stock markets where 
a -0.8 i. 
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